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ABSTRACT 

We present Very Long Baseline Array observations of the radio galaxy 3C 120 at 5, 8, 12, and 15 GHz de- 
signed to study a peculiar stationary jet feature (hereafter C80) located ~ 80 mas from the core, which was 
previously shown to display a brightness temperature ~ 600 times lager than expected at such distances. The 
high sensitivity of the images - obtained between December 2009 and June 2010 - has revealed that C80 cor- 
responds to the eastern flux density peak of an arc of emission (hereafter A80), downstream of which extends 
a large (~ 20 mas in size) bubble-like structure that resembles an inverted bow shock. The linearly polarized 
emission closely follows that of the total intensity in A80, with the electric vector position angle distributed 
nearly perpendicular to the arc-shaped structure. Despite the stationary nature of C80/A80, superluminal com- 
ponents with speeds up to 3 ± 1 c have been detected downstream from its position, resembling the behavior 
observed in the HST-1 emission complex in M87. The total and polarized emission of the C80/A80 structure, its 
lack of motion, and brightness temperature excess are best reproduced by a model based on synchrotron emis- 
sion from a conical shock with cone opening angle r] = 10°, jet viewing angle 6 = 16°, a completely tangled 
upstream magnetic field, and upstream Lorentz factor ju = 8.4. The good agreement between our observations 
and numerical modeling leads us to conclude that the peculiar feature associated with C80/A80 corresponds to 
a conical recollimation shock in the jet of 3C 120 located at a de-projected distance of ~190pc downstream 
from the nucleus. 

Subject headings: galaxies: active - galaxies: individual (3C 120) - galaxies: jets - polarization - radio 
continuum: galaxies 



1. INTRODUCTION 

The radio galaxy 3C 120 is a very active and powerful emit- 
ter of radiation at all observed w avebands. While classified 
as a Seyfert 1 ( |Burbidge]|1967| ), its spectrum shows broad 
emission lines and a complex optical morpholog y which 
has been interpreted as the result of a past merger ( [Garcia- 
Lorenzo et al.||2005] ). Variability in the emission lines has 
allowed rever beration m apping to yield a black hole mass of 
-- 3 X 1O^M0 ( Wandel et al. 1999 ). In X-rays, it is the bright- 
est broad-line radio galaxy, with the X -ray spectrum becom - 
ing softer when the intensity increases ( [Maraschi et al.|199T] ), 
suggesting that the soft X-rays are doniinated by the disk in- 
stead of the beamed jet. This has made it possible to estab- 
lish a clear connection between the accretion disk and the 
radio jet through coordinated X-ray and radio observations 
( [Marscher et al.|[2002[ [Chatterjee et al.| |2009). Along with 
3C 111, 3C 120 is one of only two broad-line radio galaxies 
detected by Fermi-LAT in the GeV photon energy range, such 
emission being most probably the result of beamed radiation 
from the r elativistic jet ob served at intermediate viewing an- 
gles (Kat aokaetal.|2011| ). 

At radio wavelengths, 3 C 120 has a blazar-like one-sided 
superluminal radio jet extending up to hundr eds of kiloparsecs 



very rich inner jet structure containing multiple superluminal 
components with apparent velocities up to 6c (Gomez e t al. 
1998, 1999i |Homan et al.|2001t[Jorstad et a l. 2005 ; MarscEer 
|et al.||2007j ). Continued monthly monitoring with the VLB A 
at 22 and 43 GHz has revealed rapid changes in the total and 
linearly polarized intensity, interpreted as resulting from the 
interaction of the jet comp onents with the external medium 



( [Gomez et al. 2060[ |2001| ). Longer wavelength Very Long 



superluminal radio jet extending up to nundreds or Kiloparsecs 
dWalker et al.fT987l [T9881 [Muxlow & Wilkinson|19 91 ). Due 
to its proximity (z=0.03 3), and to its bright millimeter emis- 
sion fAgu do et al.|2 010). Very Long Baseline Array (VLB A) 
observations at high frequencies (22, 43, and 86 GHz) show a 



^ Institute de Astrofisica de Andalucia, CSIC, Apartado 3004, 18080, 
Granada, Spain; jlgomez@iaa.es 

^ Institute for Astrophysical Research, Boston University, 725 Com- 
monwealth Avenue, Boston, MA 02215, USA 

^ School of Computing, Engineering and Physical Science, University 
of Central Lancashire, Preston PRl 2HE, UK 



Baseline Interferometry (VLBI) observations have shown ev 
idence for the existence of an underlying helical jet structure 
([Walker et al. 2001 ), whi ch has b een interpreted as the result 
oTjet precession (Hardee et al.|[2005^ Caproni & Abraham] 
[2004] ) 

Combination of the information from a sequence of 12 
monthly polarimetric VLB A observations of 3C 120 at 15, 
22 and 43 GHz has allowed imaging of the linearly polarized 
emission within the innermost ~ 10 mas jet structure, reveal- 
ing a coherent in tim e Faraday screen and polarization angles 
( [Gomez et al.|2008| ). Gradients in Faraday rotation and degree 
of polarization across and along the jet are observed, together 
with a localized region of high (~ 6000 rad m~^) Faraday 
rotation measure (RM) between approximately 3 and 4 mas 
from the core. The existence of this localized region of high 
RM, together with the observation of uncorrelated changes in 
the RM screen and RM-corrected polarization angles, suggest 
that a significant fraction of the Faraday RM found in 3C 120 
originates in foreground clouds, rather than in a sheath in ti- 
mately associated with the emitting jet (Gomez et al.' 201l| ). 

VLB A observations of 3C 120 in November 2007 revealed 
a component (hereafter C80) located at 80 mas (which cor- 
responds to a de-projected distance of ~190pc for a viewing 
angle of 16°) with a brightness temperature 7]^ ^ 5 x 10^ K, 
which is about 600 times higher than expected at such dis- 
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tances from the core ( |Roca-Sogorb et all2010| ). Analysis of 
previous observations (starting in 1982) snows that this com- 
ponent was not detected at frequencies higher than 5 GHz 
before April 2007, when was first observed at 15 GHz. Af- 
ter this epoch C80 appears in all images (even at the highest 
frequencies) at the same location withou t significant changes 
in its flux density. jRoca-Sogorb et al.] p010) conclude that 
the unusually high 7]^ of C80 could be explained by a heli- 
cal shocked jet model - and perhaps some flow acceleration 
- but it seems very unlikely that this corresponds to the usual 
shock that emerges from the core and travels downstream to 
the location of C80, requiring some other intrinsic process ca- 
pable of providing a local amplification in the density of high 
energy particles and/or magnetic field. 

Here we present the results from a new observing program 
designed to study the nature of the C80 feature. Such study is 
complemented with semi-dynamical and semi-analytical sim- 
ulations to help us interpret the structure of C80 and its sur- 
rounding jet region. The paper is organized as follows: in 
Sect. [2] we describe the new VLBA observations presented 
here and their data reduction procedures; in Sect. [3] we report 
on the results that are obtained directly from the observations; 
we present the numerical models used to reproduce these ob- 
servations in Sect. |4] Finally, in Sect. [5] we summarize our 
main results and draw the conclusions from this work. 

At the redshift of 3C 120, and under the standard ACDM 
cosmology (with Hq=11 km s~^ Mpc~^ Qm = 0.27, and Q\ = 
0.73) that we assume in this work, an angular separation of 
1 mas in the sky corresponds to a projected linear distance 
of 0.65 pc, and a proper motion of 1 mas/yr is translated to a 
superluminal speed 2.19 times larger than the speed of light. 

2. OBSERVATIONS 

The new multi-frequency and polarimetric VLBA images 
of the jet in 3C 120 presented in this paper were obtained from 
observations performed on 2009 December 14, 2010 March 
14, and 2010 June 21. The corresponding observations, that 
employed all ten antennas of the VLBA, were performed at 5, 
8, and 12 GHz (see Figs. [T}|3]), and used a 2-bit signal sam- 
pling to record in a 32 MHz bandwidth per circular polariza- 
tion at a recording rate of 256 Mbits/s. 

Calibration of the data was performed within the AIPS 
software package followin g the standard procedure for VLB I 
polarimetric observations ( Leppanen et al.|p"995| ). To cor- 
rect for opacity effects on the higher frequency observations 
(> 8 GHz) we used the recorded variation of the system tem- 
perature on every station to solve for the receiver temperature 
and zenith opacity. The final images were obtained following 
the standard hybrid-mapping procedure through an iterative 
process that emp loyed both AIPS and the DIFMAP package 
( |Shepherd|1997l ). 

For the calibration of the absolute orientation of the elec- 
tric vector position angle (EVPA) we employed both, a set of 
VLA observations of our polarization calibrators (NRAO150, 
OJ287, 3C 279, and 3C 454.3) from the VLAA^LBA Polariza- 
tion Calibration Progranj^ and the information extracted by 
the comparison of the D -terms at different observing epochs 
(see Gomez et al.||2002| ). The useful input from the first of 
these two methods was rather limited in our case because the 
VLAA^LBA Polarization Calibration Program does not in- 
clude data at 12 GHz, and the program was halted at the end 



of 2009 to focus on EVLA commissioning, hence prevent- 
ing us from getting the calibration at 5 and 8 GHz for the lat- 
est two observing epochs in 2010. However, we verified that 
the 5 and 8 GHz D-terms were very stable during our three 
observations. [Gomez et al.| ([2002 ) showed that in this case 
the D-terms can be used to calibrate those observing epochs 
without an independent EVPA calibration. This allowed us 
to calibrate the 5 and 8 GHz EVPA of the latest two observ- 
ing epochs by rotating their corresponding D-terms to match 
those of the fully-calibrated first observing epoch. 

At 12 GHz there are no available independent polarization 
data from our sources to estimate the adequate EVPA calibra- 
tion. However, we used the publicly-available total-intensity 
and linear-polarization 15 GHz VLBA images of 3C 120 
taken (and calibrated) by the MOJAVE tean]^on dates close to 
those of our multi-frequency observations (see Fig.|4]). These 
15 GHz EVPA images allowed us to calibrate our 12 GHz 
EVPA on nearby epochs by fitting Faraday rotation profiles 
at 5, 8 and 15 GHz to find the correct EVPA calibration at 
12 GHz. The stability of the D-terms at 12 GHz was also used 
as a cross-check for the final calibration at this frequency. The 
final errors in our absolute calibration of the EVPA at all our 
observing frequencies are estimated to be ~ 5°. 

3. RESULTS 

3.1. The Inner ~ 20 mas of the Jet 

Figures[T]to|4|show the new 5, 8, 12, and 15 GHz VLBA im- 
ages of 3C 120 in total and linearly polarized intensity, as well 
as their corresponding EVPA distribution. These new images 
show a straight jet within the inner r^ 20 mas that contains a 
bright and highly polarized knot at ~ 16 mas from the core 
(CI 6). The EVPA of C16 is distributed perpendicular to the 
jet axis. The same EVPA orientation is found in the C 16 coun- 
terparts at 8, 12, and 15 GHz, see Figs.[2]|4j where the images 
show an increasingly rich structure of jet knots with observ- 
ing frequency all along the innermost 20 mas of the jet. The 
multi-frequency polarization p roperties of C16 matc h those of 
the C 12 jet feature reported by |G6mez et al.| (l2011l). This, to- 
gether with the total flux dominance of both C16 and C12 in 
regions > 10 mas from the core and the fact that C16 lies in 
the jet region expected from the position of C12 and the typ- 
ical proper motions in the jet of 3 C 120 ('^ 2 mas/yr, |G6mez| 



et al.pOOT] ), allow us to identify both components as repre- 



sentmg the same jet moving knot. 

3.2. Bright Emission Region Between 80 and 100 mas from 
the Core: B80-100 

The jet feature at ~ 80 mas from the core (C80) is detected 
in all new images presented in this work from 5 to 15 GHz 
(Figs. [ip]). No other jet feature is detected in the jet down- 
stream of C80 in the 12 and 15 GHz images presented here. 
However, the new 5 GHz images show an emission region in 
an arc from north to south around the eastern side of C80, hav- 
ing the shape of an inverted bow shock, similar to those ob- 
serve d at the end of the large scale jets of FR-II radio galaxies 
(e.g., |Perley et al.| 1984b . Perhaps because of the lack of sensi- 
tivity, tETs~arc~oremission (that we have labeled as A80) was 
not detected by our previous short -integration VLBA obser- 



vations of the source at 5 GHz (see lRoca-Sogorb et al.|2010" 
Gomez et al.||2011 l, where only the detection of C80, the 



brightest region in A80, was reported 



http : //www. via .nrao . edu/astro/calib/polar/ 
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Fig. 1.— Sequence of VLBA images of 3C 120 at 5 GHz taken on 2009 December 14, 2010 March 14, and 2010 June 21 (from top to bottom). Vertical 
separations are proportional to the time difference between epochs. Total intensity contours are overlaid at 0.03, 0.07, 0.17, 0.41, 1.02, 2.49, 6.11, 14.98. 36.72, 
and 90 % of the total intensity peak at 1 .53 Jy/beam. A common convolving beam of FWHM 3.5x1.7 mas at 0° was used for all images and is shown in the lower 
left corner. Grey scale images (on a linear scale shown at the top) show the linearly polarized intensity. Black sticks (or unit length) indicate the polarization 
electric vector position angle, uncorrected for Faraday rotation. The vertical lines at ~ 80 mas from the core indicate the location where the transverse cuts of the 
fractional polarization of the jet shown in Fig.[5]were performed. (A color version of this figure is available in the online journal.) 



Together with C80/A80, the new 5 GHz images show a 
bubble-like extended emission-region larger than ~ 20 mas 
along the jet axis - from the sharp edge near the location of 
C80-, and ~ 20 mas across the jet axis. This emission struc- 
ture will be called B 80- 100 hereafter (see Fig.[T]). Neither the 
total intensity nor the polarized emission is symmetric with 
regard to the jet axis in B 80- 100. The maxima in total and 
linearly polarized intensity are both located on the southern 
side of the bubble, which resemble s the behavior of the jet in 
3C 120 within the first - 10 mas ( |G6mez et al.1[200Q| [200T] 
|2008). 

Further downstream from C80 (but still wit hin B80-100), 
C90, a mor e ext ended jet regi on reported by |Roca-Sogorb| 
et al.J ( [20TQ| ) and |G6mez et al.| p011) at observing frequen- 
cies below 8 GHz, starts to be resolved in our new images 
(Figs. [1] and [2]). Our 8 GHz images in 2009 and 2010 show 
the jet emission around C90 as a slightly more extended and 
spars er jet region than in the image taken in 2007 November 
30 by Gomez et al. (201 1 ) at the same frequency. C90 is better 
observed in our 5 GHz images, that map the entire B 80- 100 at 
higher sensitivity. This provides a reliable superluminal speed 
measurement of C90 (vc9o = 3.4 ±1.0 c), which contrasts w ith 
the stationary character of C80 ( Roca-Sogorb et al.|2010| ) lo- 
cated ~ 10 mas upstream in the jet. 



Our new images also show a bright and compact jet region 
located ~ 99 mas from the core - reported here for the first 
time - that we have labeled C99. Our kinematic study of this 
feature also reveals superluminal proper motion (vc99 = 3.0 ± 
1.1c) of similar magnitude and direction downstream the jet 
than C90. The superluminal character of the flow downstream 
from C80 resembles the behavior near the HST-1 knot in M87 
as reported by Cheung et al._ (2007j , who suggested HST-1 as 
the site of a reco llimatio n shock in the jet of M87 (see also 
IStawarz et aLp006HAsada & Nakamura|2012| ). 

Given the typical superluminal mot ions previously dete cted 
in 3C 120 on smaller jet scales (~ 4 c; 'Gomez et al.'2001), we 
consider it unlikely that C90 - as identified by Roca-Sogorb] 
[eTal. (2010) and Gomez et al. (2011 ) - can be identified with 
C99 at earlier epochs, given the exceedingly large proper mo- 
tion that would be needed (> 9c, which has not previously 
been observed in the jet of 3C 120). Our 5 GHz images show 
that C99's flux density increases progressively from 2009 De- 
cember 14, to 2010 June 21, which provides a better expla- 
nation for why C99 was not detected in our previous 5 GHz 
observations in 2007 November 30. 

3.3. Arc of Emission at ~ HO mas from the Core: A80 

As it is shown in Fig. [T] (see also Figs. [2] to |4]), the linearly 
polarized emission in A80 closely follows that of the total in- 
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Fig. 2.— Same as Fig.[T]but for the 8 GHz images. Contours levels are drawn at 0.11, 0.23, 0.49, 1.03, 2.17, 4.58, 9.64. 20.30, 42.74 and 90% of the total 
intensity peak at 1 .04 Jy/beam. A common convolving beam of FWHM 1 .8 x 0.9 mas at 0° was used for all images. (A color version of this figure is available in 
the online journal.) 
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Fig. 3. — Same as ] 



; Fig.[T]but for the 12GHz images. Contours levels are drawn at 0.17, 0.34, 0.69, 1.38, 2.76, 5.54, 11.13. 22.34, 44.83 and 90 
intensity peak at 0.75 Jy/bcc 



the online journal.) 



of the total 
eam. A common convolving beam of FWHM 1 .5 x 0.7 mas at 0° was used for all images. (A color version of this figure is available in 
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3C120 at 15 GHz 
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Fig. 4. — Same as Fig. ^ but for the 15 GHz images obtained and calibrated by the MOJAVE team (see the text). From top to bottom, the three images 
correspond to VLBA observations performed on 2009 December 10, 2010 March 10, and 2010 July 12. Contours levels are drawn at 0.13, 0.27, 0.56, 1.15, 
2.38, 4.92, 10. 17. 21 .04, 43.52 and 90 % of the total intensity peak at 0.73 Jy/beam. A common convolving beam of FWHM 1 .2 x 0.5 mas at 0° was used for all 
images. (A color version of this figure is available in the online journal.) 
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Fig. 5. — 5 GHz fractional polarization across the A80 jet structure through 
the slices (from south to north) marked in Fig.fll Diamonds correspond to 
epoch 2009 December 14, triangles for epoch 201 March 14, and squares 
for 2010 June 21. 

tensity, with values of the degree of polarization up to ~ 30% 
(see also Fig. [5]), and the EVPAs distributed nearly perpen- 
dicular to the arc-shaped structure of A80, as expected in the 
case of a shock front. 

The observations presented here show a nearly constant flux 
density evolution of C80/A80, and confirm the large bright- 
ness temperature as well as its stationarity, both within our 
new observing epochs and with regard to the reported position 
of C80 at earlier epochs. These properties appear remarkably 
similar to those expected in the case of conical recollimation 
shocks (e.g., Cawthorne & Cobb 1990; Gomez et al. 1995, 
[l997,.Agudoetal.,2001, Cawtho rne 2006 ; Nalewajko,2009 ), 
and will be discussed in detail in the next session. 

Interestingly, the region of B 80- 100 close to the jet axis 
and downstream from the sharp edge of A80 appears unpo- 
larized in our 5 GHz images. In contrast, weak linear polar- 
ization emission (with EVPA nearly perpendicular to the local 
jet axis) starts to be discerned further downstream at ~ 10 mas 
from the eastern edge of A80. The strong linear polarization 
emission of C99 shows its EVPA distribution parallel to the 
jet axis, which, together with the measured superluminal mo- 
tion, may be an indication of a travelling plane shock wave 
perpendicular to the axis. 

3.4. Faraday Rotation Screen Between 5 and 15 GHz 

Figure |6] shows the rotation measure image of 3C 120 for 
December 2009 made by combining the data between 5 and 
15 GHz. The RM images for the other two epochs of our ob- 
servations show very similar results as those for December 
2009, and are therefore not shown. As can be seen in the inset 
panels of Fig. [6J we obtain good fits to a A^ law of the EV- 
PAs throughout the jet, except at longer wavelengths in the 
innermost ~4 mas, probably due to beam depolarization and 
opacity. 

The RM screen shown in Fig. [6] is in excellent agreement 
with that obta ined for observation s taken in November 2007 
(see Fig. 7 in Gomez et al.pOll] ), providing further confir- 
mation for the stationa rity of the Faraday screen in 3C 120 
claimed previously (G omez et al.|2008 ). Note, however, that 
relatively rapid changes (in scales of tens of months) in the 
R M have also been obs erved in the innermost ~2 mas of the 
jet ( [Gomez et al.|2011 ). Our new observations also confirm 
the negative graclient m RM with distance along the jet, in- 
cluding a sign reversal at '^8 mas from the core. The motion 



of component C16 along the jet durin g the two years time 
span between the 2007 observations of Gomez et al.| ( [2011 ) 
and ours have allowed the mapping of the RM further down- 
stream, revealing a RM of -79 ± 2 rad m~^, fully consistent 
with the RM gradient along the jet mentioned previously. 

Component C80/A80 has a small rotation measure of 20 ± 2 
rad m~^, leading to a Faraday rotation in the EVPAs at our 
longest observing wavelength of 6 cm (5 GHz) of 4 degrees, 
within the estimated error in our absolute calibration of the 
EVPAs. We therefore can conclude that our EVPA maps of 
C80/A80 are not affected by Faraday rotation. 

4. A CONICAL SHOCK MODEL FOR A80. 

The structure and polarization of the A80 feature point to- 
wards an origin in the kind of shock that can occur when the 
internal pressure of a jet adjusts as a result of changes to its 
environment. This type of structure has been suggested as 
the origin of the HST-1 knot in M87 (Stawarz et al. 2006} 
,Asada & Nakamura„2012 ). A jet that is initially freely ex- 
panding may become underpressured and is therefore recolli- 
mated by a conical shock converging away from the nucleus. 
Such shocks may reflect from the jet axis, producing a second 
conical shock diverging away from the nucleus, that further 
adjusts the collimation of the jet. The result is a structure 
consisting of a pair of point-to-point conical shocks that have 
been found t o occur both in exp eriments and numerical simu- 
lations (e.g., |G6mez et al.|1995| ). In A80, there is no evidence 
for a converging shock -the shape of the structure, particu- 
larly in polarized intensity, suggests a diverging cone that will 
divert the flow away from the axis causing a drop in pressure. 
This therefore appears to indicate that A80 has arisen because 
the jet has rapidly become overpressured, perhaps as a result 
of a sudden drop in external pressure at this point. Though 
this seems the most probable explanation it is possible that 
the jet has encountered small object - perhaps a molecular 
cloud typical of the narrow-line emission regions in AGN - 
near the axis of the jet, and that the conical shock represents 
the resulting wake. This explanation seems less probable be- 
cause it depends on the lucky chance that the obstacle is near 
the jet axis, but it is in agreement with previous indications 
for a jet/clo ud collision at distances closer to the core (Gomezl 
let al.i2000| ). 

Here, the A80 feature is modelled a s a diverging conical 
shock usin g the model s desc ribed by |Cawthome & Cob"b 
( |199Qj ) and [Cawthorne] ( [2006 ). These models assume the 
plasma to be a relativistic gas with sound speed c/a/3. 

Previous monitor i ng of t he proper motions in 3C 120 (e.g., 
Gomez et al.||20001 [JOOT] I Walker et al.l[2QQTl |Homan et al.| 



2001t|Marscher et al.|2002t|Jorstad elaL 



2005] ) shows that the 



maximum apparent velocity observed is ~ 6 c. From this it is 
possible to estimate the maximum viewing angle, Omax = 19° 
and the minimum Lorentz factor of the jet, ^u = 6. (e.g., 
Gomez et al.|2000l. The value Papp = 6 is therefore assumed. 



being also relevant to the components labelled L5 to L8 by 
Walker et al. ( 2001| ) in the region upstream from C80. This 



value is used to constrain the properties of the upstream flow. 
The 5 GHz images presented in this paper are used to con- 
strain the properties of the shock wave. 

The projected semi-opening angle r]p of the conical shock is 
best estimated from the polarized intensity images presented 
earlier in this paper (Fig.[T]). The cone axis is taken to be a line 
inclined at approximately 7° from EW, pointing directly back 
to the radio core structure. The opening angle deduced from 
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Fig. 6. — Rotation measure image of 3C 120 between 5 and 15 GHz for observations taken in 2009 December 14. Contours show the 5 GHz total intensity 
image (see Fig.[l}. Black sticks indicate the RM-corrected EVR\s. Inset panels show sample fits to a A^ law of the EVPAs at some particular locations. The 
convolving beam is shown in the lower left corner. (A color version of this figure is available in the online journal.) 



the distribution of polarized flux density is about 90°, but after 
convolution of model structures with the elliptical beam, the 
best fit to this is found from a projected opening angle in the 
region of 80° corresponding to rjp = 40° . The angle between 
the polarization rods on either side of the structure (measured 
where the polarized flux density is brightest) is in the region 
of 90° . Since the model polarization E rods are orthogonal to 
the projected edge of the conical structure, a value of r]p = 40° 
corresponds to an angle of 100° between the polarization rods 
on the two sides. However, after convolution, this angle is re- 
duced to approximately 90°, in agreement with the images, 
as a result of the influence of the polarized structure nearer 
the axis. A value of r]p much less than 40° would be hard to 
reconcile with both a symmetrical cone structure and an axis 
pointing back to the core region; if the axis is maintained at 
its inclination of 7° to EW, then the southern edge of the cone 
would then exclude much of the observed flux density. Values 
of r]p significantly greater than 40° result in simulated polar- 
ization images that are inconsistent with those observed in that 
they do not reproduce the observed saddle point between the 
two sides of the polarization image. Hence, a value for r]p in 
the region of 40° seems to offer the best chance of modelling 
the observed structure of A80, and is therefore used in what 
follows. 
The true semi-opening angle r] is related to r]p by 



/3„7„sin7^> 1/V2 



(2) 



tan 7^^: 



tan 7^ 



cos I9(l-tan2 6>tan2 7^)1/2 



(1) 



which can be solved numerically for rj in terms of rjp and 0. 

A further constraint is provided by the requirement that, in 
the rest frame of the upstream flow, the shock wave must ad- 
vance toward the stationary plasma at speed greater than the 
sound speed, c/VS. This requires that 



For each value of ^, a unique model can be com puted us- 
ing the procedures described in "CawthorneJ ( |2006| ) provided 
inequality [2] is satisfied. Parameters appropriate to several 
models are given in Tablefljand the corresponding model im- 
ages have been computecTassuming that the upstream mag- 
netic field is completely disordered. The plausible range of 
is quite restricted: values below about 10° are ruled out about 
by Eq. Ej while values above 18.9° cannot yield apparent su- 
perlummal speeds as high as ~ 6c. 

The simulated images corresponding to models 1-4 resem- 
ble the observations of A80 shown in Fig. [T] The total in- 
tensity (/) images in Fig|7] show a bow wave pattern while 
the simulated polarized intensity (P) images show a similar 
pattern that is divided by a saddle point at the position of 
maximum /. The polarization E rods in the two elongated 
P features are approximately perpendicular to the direction of 
elongation, which corresponds to the outline of the conical 
shock wave. Models 1 to 4 produce simulations that differ 
mainly in two ways: first, the structures become more com- 
pact as 6 increases, because the emission from the far side 
of the shock wave becomes relatively weaker; and second, 
the degree of polarization increases as increases, mainly be- 
cause the shock becomes stronger (the compression factor of 
the shock k, decreases) as a result of increasing values of r] 
and the upstream flow velocity (^u)- 

The model best matching the degrees of polarization found 
in the images is model 3. In this model, the degree of polar- 
ization rises to about 30% near the edges of A80 in the region 
of the peak in total intensity, and levels as high as this are 
seen in this region on the southern side of A80 (as shown in 
the profiles plot of Fig.|5]). 

The simulated images for model 3 are shown in Fig. [7] The 
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TABLE 1 
Model parameters for different values of ( 



Model 




o 


Pu 


o 


hi 


o 


Pd 


Su 


^d,max 


^d.maxl^u 


Rx 


Ri 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


1 


12 


0.9874 


7.5 


0.76 


1.2 


0.9849 


4.63 


4.73 


1.02 


2.9 


3.0 


2 


14 


0.9895 


8.8 


0.50 


3.2 


0.9826 


3.62 


5.33 


1.47 


40.3 


27.5 


3 


16 


0.9929 


10.0 


0.31 


5.0 


0.9815 


2.61 


5.26 


2.02 


600 


395 


4 


18 


0.9974 


11.2 


0.14 


6.9 


0.9810 


1.40 


5.19 


3.71 


68,900 


39,300 



Note. — For all these models the upstream magnetic field was assumed to be completely tangled. 
Columns are as follows: (1) model number, (2) angle d between the axis and the line of sight, (3) /?„ upstream flow 
velocity as a fraction of c, (4) 77, the (deprojected) semi-opening angle, (5) k, the compression factor at the shock front, 
(6) ^, the angle through which the flow is deflected, and (7) /3^, the downstream flow velocity expressed as a fraction of 
c. (8) bu Doppler factor upstream the conical shock, (9) maximum Doppler factor downstream the conical shock 6d,max, 
(10) Sd,max/Su ratio of maximum downstream to upstream Doppler factor, (11) intensity ratio between the Ri downstream 
and upstream emission, and (12) a new R2 estimate of the downstream to upstream intensity with downstream intensity 
weighted by the Doppler shift. 



model is reasonably successful in reproducing the observed 
structure in the vicinity of the / peak. However, further down- 
stream the observed polarization rods turn to become orthog- 
onal to the axis whereas the model polarization rods turn in 
the opposite sense. On axis, the model polarization rods are 
parallel to the axis, but where on-axis polarization is seen in 
the observations, at distances greater than about 10 mas from 
the working surface, the rods are either orthogonal or oblique 
to the axis. It seems that if the conical shock model is ap- 
propriate, then it is the dominant influence on the polarized 
emission only within about 10 mas of the apex. 

Another way in which the model is successful is in its pre- 
diction of a very large ratio of upstream to d ownstream in- 
tensity. The observations of [Roca-Sogorb et al. (2010) have 
shown that C80 (the brightest hot spot in A80) represents an 
increase in the intensity of emission with respect to the under- 
lying jet by a factor of approximately 600, a figure that they 
found difficult to explain without invoking acceleration of the 
jet. If this is due to the formation of a shock structure, then the 
ratio of the up to downstream intensity should be about 600 
or greater. If the Doppler factors of the up and downstream 
flows are Su and S^, respectively, and the compression ratio is 
K,, then the intensity ratio should be approximately 



R 



= ^A = ^(5a-6)/3 



Sd 



(2-a) 



where a ::^ -1 is the spectral index for C80 ( |Roca-Sogorb et al. 
2010| ). In estimating the value of Ri the largest downstream 
Doppler factor S^^max (for flow at angle 0-^ to the line of 
sight) has been assumed. The values of k,, Sd,max for the down- 
stream flow closest to the line of sight and Ri are shown in Ta- 
ble [T] Because 6^ < S^^max for most of the downstream flow, 
Ri IS likely to be an overestimate of the intensity ratio, and so 
a second estimate based on an intensity weighted average of 
the downstream Doppler factor, 7?2, has also been included in 
the table. 

Table [T] shows that, as increases from 12° to 18°, so the 
downstream to upstream brightness ratios, Ri and R2 increase 
very rapidly. The greater part of this effect is due to the com- 
pression of the plasma (indicated by the value of hz), and its 
effect on the magnetic field and particle density. However 
the increase in the value of Sd^max/^u is also significant. This 
occurs because (i) Su is decreasing because is increasing 
and (3 is also increasing in the regime where sin > %^ ; and 
(ii) because S^ is varying only weakly, which occurs because 
although f3d decreases and Od increases, the shock becomes 



stronger and parts of the downstream flow are directed closer 
to the observer's line of sight. 

From Table [T| it is clear that parameters similar to those of 
model 3 (6 = lo°), which gave the closest match for the frac- 
tional polarization values, should have a downstream to up- 
stream brightness ratio of order several hundred. It therefore 
seems likely that such models can reasonably reproduce both 
the observed polarization structures and the very high rati o of 
down to upstream intensity demanded by Roca-Sogorb et al.| 
( |20TQl ). 

A few variants on the models described here have been in- 
vestigated to determine whether they might explain the trans- 
verse polarization rods seen on the downstream side of C80. 
First a poloi dal magne t ic fiel d component was added, as de- 
scribed by Cawthome| p006| ). The influence of the poloidal 
field on the polarization is strongest on the near side of the 
shock front, producing transverse polarization rods near the 
peak in /, where none are seen. The effect is negligible on 
the far side of the shock, which extends further downstream 
toward the region where transverse polarization rods would 
be desirable. The possibility that the upstream flow might be 
converging slightly has also been investigated, as suggested 
by ,Nalewajko, ( ,2009^ . The effect was to increase the degree 
of polarization at a given value of 6, but otherwise the sim- 
ulated images were very similar. It seems that neither of 
these modifications helped to explain the larger scale polar- 
ization properties beyond A80. It is possible that allowing for 
a more complex upstream field structure might explain the en- 
tire polarization structure of B 80- 100, but on the other hand, 
it seems likely that the conical shock is not the dominant in- 
fluence on polarization toward the downstream edge of the 
B 80- 100 structure. Indeed, the superluminal motion found 
for C90 and C99 suggests that instead the polarization struc- 
ture downstream of A80 may be the result of moving plane- 
perpendicular shock waves that appear on the wake of the A80 
conical shock. 

It is worth noting that, although the upstream Doppler factor 
is rather low in model 3 (::^ 2.6) the jet to counter-jet ratio is 
still reasonably high, giving 



^counter jet 



l+PuCO^O 

l-/3u COS 6 



2-a 



1.2 xlO"^ 



for a !::^ -0.5, as expected nearer the core. 



5. SUMMARY AND CONCLUSIONS 



A Recollimation Shock 80 mas from the Core in the Jet of the Radio Galaxy 3C 120 
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Fig. 7. — Simulated images and plots corresponding to model 3 (see Table^. Plots of the six panels are as follows: A) outline of the shock structure and the 
(unconvolved) polarization rods of length proportional to P , B) contours of/ aner convolution with a beam chosen to match in shape that of the observations; the 
convolved polarization is shown by rods parallel to the E field and of length proportional to P, C) contours of convolved P, with polarization angle x shown by 
the orientation of rods of constant length, D) variation of fractional polarization (from the convolved images), E) along the jet axis, and F) across the jet through 
the region of the / peak. In all six plots the units of distance are scaled so that the height of the cone (from apex to base) is unity. 



First anal ysis of the properties of com ponent C80 were car- 
ried out by |Roca-Sogorb et aL ( |2010| ), concluding that al- 
though a helical shocked jet model - including perhaps some 
bulk flow acceleration - could explain the unusually large 
brightness temperature, it appears unlikely that it corresponds 
to the usual shock that moves from the core to the location of 



C80. Rather, Roca- Sogorb et al. (2010) proposed the need for 
an alternative process capable of explaining the high bright- 
ness temperature of C80, its appearance in high frequency im- 
ages after April 2007, and its apparent stationarity. One of the 
possible alternatives considered by these authors involved a 
stationary shock at the location of C80 with emission sud- 
denly enhanced by the arrival of a region of enhanced particle 
number in the jet flow, such as the one proposed to explain the 
kinematic and flux evolut ion properties of the HST-1 knot in 
M87 (Sta warzetal.|2006| ). 

The main goal of the new observations presented here is 
to further constrain the physical processes taking place in the 
jet of 3C 120 that have led to the extreme properties of C80. 
Thanks to the increased sensitivity achieved by our new ob- 
servations we report on the existence of the A80 arc of to- 
tal intensity and linear polarization associated to C80, as part 
of a larger bulge of emission that extends ~20 mas along 



and across the jet. Most importantly, the polarization vec- 
tors are distributed perpendicularly to the semi-circular shape 
of A80, as it would be expected for the case of a compres- 
sion by a shock front. This evidence, together with the ex- 
cess in brightness temperature displayed by C80/A80 and its 
stationarity in flux and position support the model suggested 
by Roca-Sogorb et al. (2010) in which C80/A80 corresponds 
to a standing shock. More specifically, our modeling sug- 
gests a conical recollimation shock. Indeed, our simulations 
based on the synchrotron emission from a conical shock, as 
described by Cawthome (2006), reproduce quite closely the 
observed total and linearly polarized emission structure, the 
electric vector distribution, and the increased brightness tem- 
perature of C80/A80, allowing constraints on the values of the 
jet flow in 3C 120 and the geometry of the conical shock at 
~ 80 mas from the core. In particular, our simulations pro- 
vide the cone opening angle r] = 10°, the jet viewing angle 
6> = 16° at the location of A80, and the upstream Lorentz fac- 
tor 7^ =8.4. 

An important issue to investigate further is the origin of the 
recollimation shock at such large distances from the core of 
emission in relativistic jets as in 3C 120, and presumably in 
M87. While for M87 it has been proposed that the recolli- 
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mation origina tes from the transition betw een a parabolic to a 
conical shape f Asada & Nakamura|2012| ), in 3C 120 we sug- 
gest that the most plausible cause is a sudden drop in the ex- 
ternal pressure, leading to the formation of a conical shock 
wave opening away from the nucleus. 

The study of the nature of C80 - and the structure of its re- 
lated A80 arc of emission - that we present in this paper were 
only possible thanks to the high angular resolution provided 
by the VLBA and its large sensitivity and good performance 
for polarimetric observations. The use of synthetic images of 
the total intensity and linear polarization of conical shocks has 
also proven to be a powerful tool to interpret the nature of jet 
structures, and to constrain the physical and geometrical prop- 
erties of such structures and the jet plasma that forms them. 
Studies similar to that presented here can be carried out for a 
number of cases - additional potential examples might be the 



HST-1 and Kl knots in M8 7 ( [Stawarz et al.|2006| ) and 3C 380 
( [Papageorgiou et al. |2006| ) . Performing these studies will be 
important to obtain relevant information on the relativistic jets 
and their surrounding medium. 
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